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HIGH- RESOLUTION RADIO OBSERVATIONS OF THE REMNANT OF SN 1987A AT HIGH FREQUENCIES 
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ABSTRACT 

We present new imaging observations of the remnant of Supernova (SN) 1987A at 44 GHz, performed 
in 2011 with the Australia Telescope Compact Array (ATCA). The 0'/35 x 0'.'23 resolution of the 
diffraction-limited image is the highest achieved to date in high-dynamic range. We also present a 
new ATCA image at 18 GHz derived from 2011 observations, which is super-resolved to 0'.'25. The 
flux density is 40±2 mJy at 44 GHz and 81±6 mJy at 18 GHz. At both frequencies, the remnant 
exhibits a ring-like emission with two prominent lobes, and an east- west brightness asymmetry that 
peaks on the eastern lobe. A central feature of fainter emission appears at 44 GHz. A comparison 
with previous ATCA observations at 18 and 36 GHz highlights higher expansion velocities of the 
remnant eastern side. The 18-44 GHz spectral index is a = —0.80 (S v oc v a ). The spectral index 
map suggests slightly steeper values at the brightest sites on the eastern lobe, whereas flatter values 
are associated with the inner regions. The remnant morphology at 44 GHz generally matches the 
structure seen with contemporaneous X-ray and Ha observations. Unlike the Ha emission, both the 
radio and X-ray emission peaks on the eastern lobe. The regions of flatter spectral index align and 
partially overlap with the optically-visible ejecta. Simple free-free absorption models suggest that 
emission from a pulsar wind nebula or a compact source inside the remnant may now be detectable 
at high frequencies, or at low frequencies if there are holes in the ionised component of the ejecta. 

Subject headings: circumstellar matter — radio continuum — supernovae: individual (SN 1987A) - 
supernova remnants — acceleration of particles — radiation mechanisms: miscel- 
laneous 



1. INTRODUCTION 

Radio supernovae (SNe) result from the collision be- 
tween a supernova shock and the progenitor's circumstel- 
lar medium (CSM). As the interaction between the prop- 
agating blast wave and the CSM drives the particle accel- 
eration process, supernova remnants (SNRs) are natural 
laboratories for studying particle sp ectra and their varia- 
tion in time (see lWeiler et al. 1120021 for a review). 

Supernova 1987A in the Large Magellanic Cloud, as 
the only nearby core-collapse supernova observed to 
date, has provided a unique opportunity to study 
the evolution of the SN-CSM interaction in young ra- 
dio remnants. The complex CSM distribution in the 
SNR is believed to have originated from a red su- 
pergiant (RSG) which has evolved into a blue su- 
pergi ant (BSG) about 20 000 y ears before the explo- 
sion (jCrotts fe Heathcote 1 120001 ). Models of the pro- 
genitor evolution suggest that the equatorially denser, 
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i.e. slower, RSG wind (|Blondin fe Lundavist. I 119931 : 
iMartin fc Arnett II1995D was swept up by the faster BSG 
wind orris fc Pod siadlowski 2007), thus forming high 
density rings. In particular, beside the central circular 
ring in the equatorial plane [equatorial ring, ER), obser- 
vations with the Hubble Space Telescope (HST) have also 
revealed two outer rings that formed from the mass loss 
of the progen itor star, located on either side of the equa - 
torial plane (|Jakobsen et al. 1 119911 iPlait et all 119951) , 
which confer to SNR 1987A a peculiar triple-ring nebula 
structure. 



Since the radi o dete ction of the remnant in mid-1990 
(|Turtle et al.. I I1990D . the synchrotron emission has 
been generated by the shock wave propagating into 
the ring-shaped distribution of the CSM in the equa- 
torial plane. Monitoring of the flux density has been 
regularly undertaken with the Molonglo Observatory 
Synthesis Telescope (MOST) at 843 MHz and at 
1.4, 2.4, 4.8 and 8.6 GHz wit h the Australia Tele- 
scope Compact Array (ATCA) (|Stavelev-Smith et al. I 
1992). A TCA observations have b e en ongoing for 
25 ye ars dStavelev-Smith et alJll993t iGaensler et al. I 
20011: Manchester et al.. I l2002t 



19971: IBall et al. 
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2002; 

iZanardo et alT 12010 ). 

An exponential increase of the flux density has been 
measured at all frequencies since day ^5000 after the 
explosion, which is likely due to an increasing efficiency 
of the acceleration process of particles by the shock front 
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(jZanardo et al. II201QT ). 

The morphology of the non-thermal radiation emitted 
by relativistic electrons accelerated in the remnant, has 
been investigated using im ages at 9 GHz since 1992 
(jStaveley-Smith et~ al. 1993), with a spatial resolution 
of Of.' 5 achieved via maximum entropy super-resolution 
(jGaensler et al. I Il997t iNg et al. I 120081 ) . These images 
have provided the first insight into the marked east-west 
asymmetry of the radio emission. As receivers operat- 
ing in the 12 mm band (16-26 GHz) were introduced 
at the ATCA in 2001, the first imaging observations at 
18 GHz were unde rtaken in 2003 July, at an effective 
resolution of 0"45 ([Manchester et al. 1 120051 ) . Very long 
baseline interferometry (VLBI) obs ervations of the SNR 
were successful in 2007 October (iTingav et al. 1 120091) 
and 2008 November (|Ng et al. 1201 ID at 1.4 and 1.7 GHz, 
respectively. These observations provided the first im- 
ages with ~ O'.'l resolution, but with low sensitivity and 
dynamic range. Nevertheless, the VLBI images captured 
the presence of small-scale structures in bright regions 
(]Ng et al. 1 1201 lh . At the same time, ATCA observa- 
tions at 36 GHz in 2008 April and October resulted in 
high -dynamic range ima ges with an angular resolution of 
0'.'3 ([Potter et al. 1 120091 ) . The 36 GHz images, in com- 
bination with the 9 GHz observations at similar epochs, 
were used to construct the first spectral index image of 
the SNR, with a resolution of 0'.'45. The resultant map 
provided the first glimpse into the spatial variations of 
the synchrotron spectral index across the remnant, and 
hinted at the possible association between flatter spec- 
tral indices and the bright sites in both lobes. After 
the ATCA upgrade in mid-2009 with the Compact Ar- 
ray Broadband Backend (CABB) (Wilson et al. 2011), 
the remnant was imaged at higher frequencies. The first 
resolved image at 94 GHz was produc ed from observa- 
tions between 2011 June and August (jLakicevic et al. I 
l2012t ). Comparison of the 94 GHz image with data at 9 
GHz yielded a low-resolution spectral index map, with a 
much larger region of flatter spectral index on the eastern 
lobe than on the western lobe. 

This paper presents the first high-resolution image of 
SNR 1987A at 44 GHz and a new image at 18 GHz, both 
derived from ATCA observations performed in 2011. In 
§ [2] and § [3] we describe the imaging procedure and the 
resultant remnant morphology. In § 21 we estimate the 
remnant expansion rate, by comparing the 44 GHz data 
with the 36 GHz data from 2008, while the new 18 GHz 
data are compared to the 2003 datasets at 17 and 19 
GHz. In §[5] we assess the 18-44 GHz spectral index dis- 
tribution, via different spectral mapping techniques. In 
§ [51 the comparison of the 44 GHz image and contempo- 
raneous X-ray and Ha observations, is discussed. In § [7] 
we assess the likelihood of detecting at 44 GHz the ra- 
diation emitted by a compact source located within the 
inner regions of the remnant. 

2. OBSERVATIONS 

SNR 1987A was observed at 44 GHz with the ATCA 
in 2011 January and November. The observations were 
performed on January 24 with the array in 6A config- 
uration and on November 19 with the array in 1.5D 
configuration, with maximum baselines of 5939 m and 



4439 m, respectively. The observations were taken over 
2x2-GHz bandwidth and centred on 43 and 45 GHz. 
The observations at 18 GHz were performed on 2011 
January 26 with the ATCA in 6A configuration, in 
two bands, each of 2-GHz width, centred on 17 and 
19 GHz. Atmospheric conditions were exceptional dur- 
ing both the 44 GHz sessions with rms of the path 
length fluctuations below 250 /zm, thus confirming that 
phase stability conditions at the ATCA can be met not 
only in winter but also in spring an d summer nights 
(Mid delberg. Sault fc Kesteven I [20061 ) . The rms of the 
path length fluctuations was below 300 /Ltm during the 
18 GHz session. The parameters of the three observing 
sessions are listed in Table [TJ 

In all observations, the standard bandpass calibrator 
PKS B0637-752 was observed for 2 minutes every 90 
minutes, while the phase calibrator PKS B0530-727 was 
observed for 1.5 minutes every 6 minutes on the source. 
Uranus was used as flux density calibrator at 43 and 45 
GHz. At 17 and 19 GHz, we used PKS B1934-638 as 
primary flux density calibrator. This is tied to Mars and 
ha s a stable flux density at 12 mm- wavelength, as shown 
in ISault I (|2003D . The absolute flux scale of both Mars 
and Uranus between 1 and 50 GHz has been shown t o 
have an uncertainty of 1 — 3% (|Perlev fc Butleri r2012'). 
With good atmospheric conditions and good calibration, 
we estimate that the flux calibration accuracy at 18 and 
44 GHz is within 5%. The observations were centred on 
RA 05 h 35 m 27!975, Dec -69° 16' IT/08 (J2000), as in 
Potter et al. (2009). 

The MIRIAeQ data reduction package was used to pro- 
cess all datasets. For the 44 GHz observations, the 
task ATFix was first used to apply corrections to the 
system temperatures, instrumental phases and baseline 
lengths. The Ja nuary and Nove mber datasets were im- 
aged, CLEAN-ed (|Hogbom Ill974[ ) and self-calibrated sep- 
arately before being combined. In particular, the prelim- 
inary CLEAN model was constructed by using 400 itera- 
tions for the January dataset and 180 for the Novem- 
ber dataset. Phase self-calibration was performed on 
the separate datasets and on the combined data, over 
a 2-minute solution interv al. For imagin g, a weighting 
parameter of robust = 0.5 (|Briggs. Ill995| ) was used. De- 
convolution of the combined data was performed with 
the m aximum entropy method (MEM) (|Gull fc Daniell I 
I1978D . The MEM model of the combined data led to 
diffraction- limited image resolution of 0'.'35 x 0'.'23 with a 
moderate dynamic range (~165). The diffraction-limited 
image was then slightly super-resolved using a 0'.'25 cir- 
cular beam. Figure Q] shows the combined Stokes-I con- 
tinuum image derived from the January and November 
observations before and after super-resolution. The im- 
age yields an integrated flux density for the SNR of 40 ± 2 
mJy at day 8886 after the explosion, where the rms error 
in the image and the flux density uncertainty are added 
in quadrature. In Stokes-Q, U and V images, the source 
is not detected. The 3<r upper limits to the flux densities 
are 60 ^tJy beam -1 or 2% of the maximum Stokes-I flux 
density. 

The data from the 18 GHz observations were flagged, 

1 http:/ /www. atnf.csiro.au/computing/software/miriad/ 
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TABLE 1 

Observing Parameters 



Parameter 


18 GHz 


44 GHz 


Date 


2011 Jan 26 


2011 Jan 24 


2011 Nov 19 


Day since explosion 


8738 


8736 


9036 


Center frequencies* (GHz) 


16.963 & 18.964 


43.026 & 45.026 


No. of antennas 


6 


6 


6 


Array configuration 


6A 


6A 


1.5D 


Total observing time(hr) 


9.18 


9.68 


9.82 


Averaged rms path length (/nn) 


274 


236 


198 



* CABB wideband mode on two interleaved 2.0-GHz wide frequency bands, each 
with 2048 X 1-MHz channels. 



then split into the separate observing bands until the 
imaging and deconvolution steps, so as to account 
for frequency-dependent terms in the calibration step. 
Phase self-calibration was performed over a 2-minute so- 
lution interval, which gave sufficient S/N ratio. As for 
the reduction procedure of the 44 GHz data, the weight- 
ing parameter robust was set at 0.5 and, after data de- 
convolution via MEM, the resulting image was restored 
with the diffraction-limited beam of half-maximum size 
0'.'63 x 0'.'47. Because of the high S/N, a super- resolved 
image was also obtained by r estoring the MEM mode l 
with a 0'.'25 circular beam, as [Manchest er et af~l ([2005D 
have shown that this degree of super-resolution gives re- 
liable images. The diffraction-limited and super-resolved 
versions of the Stokes-I continuum 18 GHz image of SNR 
1987A, as derived at day 8738 since the explosion, are 
shown in Figure [2] The image gives an integrated flux 
density for the SNR of 81 ±6 mJy. The image parameters 
are summarised in Table [5J 

3. MORPHOLOGY OF THE RADIO EMISSION 

3.1. 44 GHz morphology 

The emission at 44 GHz appears to be mainly distributed 
in an elliptical ring, with brightness peaking on the east- 
ern lobe. The east-west asymmetry of the radio emis- 
sion is a characteristic of SNR 1987A, which has em erged 
since the first ATC A images (jGaensler et al. lfl997T) , and 
has been monitored via the super-resolved ATCA images 
at 9 GHz (|Ng et al. I [20081 hereafter N08). This asym- 
metry has been confirmed in ATCA obse r vations at 18, 
36 and 94 GHz ([Manchester et al. 1 120051: IPotter et al~1 
120091 : lLakicevicet al. Il2012t ). Slices through the 44 GHz 
image at six position angles are shown in Figure [3l The 
brightness of the eastern lobe peaks between the radial 
profiles at 60° and 90°, while the ratio of the profile max- 
ima corresponding to the brightness peaks of the eastern 
and western lobes, is ~1.5. Overall, the eastern half of 
the image appears significantly brighter than the west- 
ern one, with a ratio of ~1.6 between the integrated flux 
densities of the two regions, i.e. east and west of the 
geometric centre, approximately located ^75 mas east 
of the VLBI position of the SN (jRevnolds et al. 1119951) 
[RA 05 h 35 m 27!968, Dec -69° 16' 11'.'09 (J2000)]. This 
emission ratio is higher than the ~ 1.4 obta ined at 9 GHz 
(N08 ) and the ~ 1.3 measured at 36 GHz (jPotter et al. I 
2009), via fitting the images with an equatorial belt torus 
characterised by a west-east linear gradient of the bright- 
ness distribution (see N08). 



The profiles also show a possible third peak slightly west 
of the centre of the remnant, which reaches maximum in 
the 120° slice, at 47% of the brightness peak of the west- 
ern lobe. The overall emission from this central region is 
1.4 ± 0.2 mJy, thus making ~ 4% of the total integrated 
flux density. 

In Figure [4] the 44 GHz image is compared to the 2008 
observations at 36 GHz (see Table [3]). The remnant ex- 
pansion over the 4- year time frame is noticeable from 
the contours of the peak flux density at identical levels. 
A total integrated flux density of 23.6 ± 5.2 mJy is de- 
rived by scaling the flux density of the 36 GHz image 
to that associated with the remnant at 44 GHz at the 
same epoch, as a = —0.68 results at day 7815 from the 
spectral index function that fits observations betw een 1.4 
and 8.6 GHz (see Figure 7 in lZanardo et al. Il2010t ). This 
leads to a yearly flux density increase of 24 ± 6%, which 
is higher than the 17 ± 8 % per year derived at 8.6 GHz 
bv lZanardo et alTl (|2010D . from monitoring observations 
between Year 20 and 21. 

In Figure the 44 GHz image is compared with that 
derived by fitting the N08 inclined-torus model to the 
9 GHz observations performed in 2011 April 22 (Ng et 
al., in preparation). The fitted model has a radius of 
0'.'917, which surpasses by ^67 mas that u sed to fit the 
36 GHz image in 2008 (potter et al. Il2009[ ). The flux of 
the 9 GHz model is scaled to that corresponding at 44 
GHz, via the spectral index a = —0.74, which is extrap- 
olated from contemporaneous ATCA observations span- 
ning from 1.4 to 94 GHz (see Figure [6]). The resultant 
model image, as shown in FigureEl has a total integrated 
flux density of 40 ± 1 mJy and the restoring beam of the 
diffraction- limited image at 44 GHz. The residual be- 
tween the model and the actual visibilities demonstrates 
that the model matches the observations along the pro- 
file at PA 90°, where the ratio of the brightness peaks 
between the two lobes equals the ~ 1.6 measured in the 
observations. The residual visibilities peak on the east- 
ern lobe at PA 30° and 150°, since the new observations 
exhibit more extended bright regions than what repro- 
duced by the linear gradient of the modelled flux distri- 
bution. The linear gradient assumption also constrains 
the ratio of the integrated flux density over the eastern 
and western halves of the image, to ~ 1.3. Therefore, the 
N08 model, which fits well the radio remnant as seen at 
9 GHz, might not be applicable to the large-scale emis- 
sion morphology that emerges from the 44 GHz observa- 
tions. 
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TABLE 2 
Image Parameters 



Parameter 


44 GHz 


18 GHz 


Integrated flux density (mjy) 


40 ±2 


81 ±6 


Restoring beam (") 


0.350 x 0.225 


0.630 x 0.474 


Position angle (°) 


1.674 


0.879 


Rms noise (mjy/beam) 


0.023 


0.088 


Dynamic range 


~ 165 


~ 176 



3.2. 18 GHz morphology 



In the diffraction-limited image, the emission at 18 GHz 
looks primarily ring-shaped, while a secondary struc- 
ture seems to emerge from the super-resolved image (see 
Figure [5]). Two arm- like features appear to stem from 
the northern and southern edges of the major ellipsoidal 
structure of the emission, peaking at PA 0°, and extend- 
ing towards the western side of the remnant. While these 
features suggest a double-ring formation within the SNR, 
they might more likely be noise artefacts of the super- 
resolution process. 

In Figure the super-resolved image is compared to 
that derived from previous ATCA observations at 18 
GHz performed in 2003 July 31, which w as also super- 
resolved to 0'.'25 ([Manchester et al. | |2~00& Table© . The 
morphology similarities between the 2003 and 2011 im- 
ages emphasise the significant expansion of the remnant 
over 8 years. Slices through the 2011 image indicate 
that the central emission varies from 20% to 24% of 
the brightness on the eastern lobe. This figure is consis- 
tent with the 2003 observations, where the central emis- 
sion was estimated —20% of the average ring intensity 
(jManchester et al. I [20051 ) . 

The integrated flux density of the 2003 image, as derived 
from the combination of the data at 17.3 and 19.6 GHz, 
was 27 ± 4 mjy, which is lower than what expected from 
an extrapolation of the radio spectrum. This value leads 
to a threefold increase in flux from day 6003 to 8738, 
i.e. a yearly increase rate of 27 ± 5%. This increase rate 
exceeds the 17 ± 8% per year derived at 8.6 GHz by 
iZanardo et al. I (|201Q[ ). for data between Year 20 and 21 
since the SN. 

In terms of remnant asymmetry, at 18 GHz the radial 
profiles in Figure [7] show that the asymmetry peaks 
on the eastern lobe between PA 60° and 90°. The 
asymmetry ratio between the eastern and western 
brightness peaks is ~ 1.4, and a similar ratio is obtained 
from the ratio between the flux density integrated over 
the eastern and western halves of the image. This value 
matches the asymmetry ratio derived by N08 from 
1992 - 2008 data at 9 GHz. 

4. REMNANT EXPANSION 

From 1987 to 1992, i.e. in the early stages of the SNR, 
the remnant had an expan sion velocity of ~ 30,000 km 
s _1 (|Gaensler et al. Ifl997t l. When the radio emission re- 
emerged, this was followed by a d rastic deceleration to 
3000 km s" 1 (|Gaensler et al. Ill997l) . From 1992 to 2008, 
the expansion has been determined from the radius of the 
9 GHz images fitted via a torus model, with a resultant 



velocity of 4000±400 km s" 1 from -day 1800 to day 7620 
(N08). 

In Figure [9l the expansion of the remnant at 18 and 44 
GHz is estimated as the change over ti me of the distance 
betw een the VLBI position of the SN ([Reynolds et al. I 
Il995h and the mean position of the peak of brightness 
for the eastern and western lobes, as derived from the 
emission profile at PA 90° centred on the SN coordi- 
nates. Comparison between observations at 18 GHz at 
day 6003 and 8738, gives a velocity of 4100 km s _1 (blue 
line), while the comparison between the 36 and 44 GHz 
datasets, within the 2008-2011 time frame, leads to a ve- 
locity of 3900 km s _1 (purple line). An identical result 
of 3900 ± 300 km s" 1 is obtained if the 2011 dataset at 
18 GHz is included in the fit of data between 2008 and 
2011, since the 18 GHz measurements align with those 
at 36 and 44 GHz. All fits are consistent with the 9 GHz 
results. 

With reference to the SN position, it should be noted 
that both the datasets at 18 GHz (2003-2011), and the 
datasets at 36 GHz (2008) and 44 GHz (2011), show an 
asymmetric expansion of the remnant, with a larger in- 
crease of the distances between the SN position and that 
of the brightness peak on the eastern lobe, compared to 
the western side. The asymmetry in the expansion has 
been noted since the images at 9 GHz in early epochs 
(jGaensler et al. 1120071 ) . where an asymmetry in the ini- 
tial expansion of the SN ejecta has been proposed as the 
explanation for the fact that the eastern lobe is located 
further than the western one from the SN site. In the 
case of the 18 GHz datasets, an expansion velocity of 
5900 km s" 1 is derived for the eastern lobe, while 2300 
km s" 1 is the expansion velocity of the western bright- 
ness peak. Similarly, in the datasets at 36 and 44 GHz, 
the eastern expansion velocity is 6300 km s _1 while the 
western is 1600 km s . Averaging the values, the ex- 
pansion velocity of the remnant western lobe is 1900±400 
km s _1 , while that of the eastern lobe is 6100±200 km 
s -1 . While these values are affected by some uncertainty 
in identifying the location of the remnant with respect to 
the VLBI coordinates of the SN site, the results from the 
18, 36 and 44 GHz data imply that the expansion veloc- 
ity of the eastern lobe is significantly larger than that of 
the western lobe. Therefore, a marked asymmetry of the 
remnant emerges not only in terms of brightness distri- 
bution of the radio emission but also in terms of remnant 
expansion. 

5. SPECTRAL INDEX MEASUREMENTS 

Spectral index maps allow the examination of the spa- 
tial variations of the spectral index, a (where S oc i/ Q ), 
across the remnant and, therefore, of the electron accel- 
eration processes associated with the propagation of the 
blast wave. The 2011 observations at 18 and 44 GHz 
have been used to derive a map of the distribution of 
the spectral index in the SNR. To match the u — v cov- 
erage of the 18 GHz observations, the visibilities from 
baselines greater than 400 kA were filtered out in the 
44 GHz datasets. At both frequencies, the images have 
been derived using an identical reduction procedure. In 
particular, identical weighting parameters, deconvolution 
algorithm and phase self-calibration iterations were ap- 
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TABLE 3 

Pre-CABB ATCA Images at 18 and 36 GHz 

Parameter 18 GHz a 36 GHz 6 

Observing session 2003 July 31 2008 Apr 25 & Oct 7,8,12 

Day since explosion 6003 7815 c 

Centre frequencies (GHz) 17.3 & 19.6 34.9 & 37.4 

Flux density (mjy) 27 ± 4 27 ± 6 

Restoring Beam (") 0.45 X 0.39 0.33 X 0.24 

Position Angle (°) 2° -1.3° 

Super-resolution (") 0.25 — 

a Sec Manchester ct al. (2005). 

b See lPotter et al. I RMM . 

c Average date of the combined datasets. 



plied to both datasets. It should be noted that the self- 
calibration technique, while improving the image resolu- 
tion, removes all absolute positional information. There- 
fore, the astrometry of each image was compared with 
that prior to self-calibration and the images were shifted 
to align with recent VLBI observations of the remnant 
(Zanardo et al., in preparation). This process allowed 
us to accurately place the position of the remnant with 
respect to the VLBI position of the SN mentioned in § G3 
The aligned images were then restored with a 0'/4 circular 
beam and regridded at a pixel scale of 3 mas. 

5.1. Methods for measuring the spectral index 

To verify that any claimed spectral index variations have 
not been artificially generated by the data reduction 
and/or are tied to a particular spectral mapping tech- 
nique, the spectral index distribution within the remnant 
is reconstructed via three different methods: (z) spectral 
tomography (at); (ii) flux ratio (as); (Hi) temperature- 
temperature (or T-T) plots (oitt)- 

With the method of spectral tomography 
(|Katz-Stone fc Rudnick I I1997D . a difference image, 
I t , is calculated by scaling the 44 GHz image by a trial 
spectral index, at, and subtracting it from the 18 GHz 
image, as it follows: 

^44 

where lis an d I44 are the images at frequency vis and 
^44, respectively. When a t reaches the actual spectral 
index of a particular feature, the feature vanishes in 
the local background of the difference image. On the 
other hand, if a component has spectral index greater 
or smaller than at, the difference image appears with a 
distinctively positive or negative residual compared to 
the local background, as this component will get over or 
under-subtracted. 

The second method used to construct a spectral in- 
dex map is that of direct image division, being as = 
log(5i8/<S44)/log(18/44). As this method is sensitive 
to variations in background/foreground emission that 
might cause spatial changes in the observed spectral in- 
dex, it is not usually the preferred method to inves- 
tigate the spectral index dis tribution over large SNRs 
(jAnderson fc Rudnick. 111993ft . 

To validate the results from the flux density ratio 
Sxs/Si4,, the spectral variations over small regions within 



the remnant are also assessed from the flux slope m, 
wh ere Sir = wlSaa + q- This approach was first applied 
by iCostain I (|196fl and iTurtle et al~l (|1962D to data of 
very low angular resolution, where instead of fluxes, the 
brightness temperatures, Tb oc v~( a+2 \ are plotted at 
two frequencies, and the spectral index is derived from 
the best fitting line of the plotted points (thus referred 
to as temperature-temperature or T-T plot). 

The slope of the best-fitting line to the Tf, points derived 
for each pixel of the images at 18 and 44 GHz (Tig and 
T44), within a specific region, yields the spectral index 
a.TT of that region. While small differences in the back- 
ground levels cause spreading of the plotted T-T points, 
because of the intrinsic spatial averaging of the linear 
regression method, the regional oltt values are less af- 
fected than as by zero- level data (jLeahy fc Roger lfl99lL 
iZhang et allfl997ft . 

5.2. Interpretation of spectral index variations 

Understanding the spatial variation of the spectral index 
and emission across young SNRs is important to inves- 
tigate the structure of the expanding shock. The spec- 
tral index can be used to probe the compression ratio of 
the shock on a local scale. In turn, the emission probes 
cosmic ray and magnetic field density. Combined, these 
provide the basis for investigating the nature of shock 
acceleration processes. 

In the tomography gallery in Figure [TOJ spectral in- 
dex variations are investigated in the range between 
— 1.4 < at < —0.3, with increment 5at = 0.1. It can be 
seen that the spectral index distribution across the rem- 
nant is not uniform. In detail, the inner areas of the rem- 
nant located along PA ~30° (210°), appear to match the 
grey background for —0.6 < at < —0.3, with at ~ —0.4 
in the inner part of northern region, and a t ~ —0.5 in 
the inner part of the feature located slightly SW of the 
SN position. The spectral index gradually steepens to- 
wards the edges, with the northern and southern regions 
becoming over-subtracted at at ~ —0.7 and at ~ —1.0, 
respectively. For the bright regions on the eastern lobe 
of the remnant, sites within the contours at 90% of the 
peak flux density (see Figure 1 1 0j> primarily match the 
grey background for at ~ —0.8. Within the 60% con- 
tours, a t changes from —0.7 (southern part), to —0.8 
(central region) and —0.9 (northern region). The en- 
tire bright area on the eastern lobe is over-subtracted 
for at ~ —1.1. On the western lobe, sites within the 
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TABLE 4 








Regional Spectral indices 




.Hegion 




a s„ c 


Aa d 


Aa/ag^ 










(%) 


Nl 


-0.88 ±0.13 


-0.67 ±0.05 


0.21 


31 


N2 


-0.65 ±0.25 


-0.71 ±0.05 


0.06 


8 


El 


-0.91 ±0.07 


-0.79 ± 0.06 


0.12 


16 


E2 


-0.80 ±0.06 


-0.66 ±0.05 


0.14 


22 


SI 


-0.85 ±0.09 


-0.73 ± 0.05 


0.12 


16 


S2 


-0.72 ± 0.09 


-0.85 ± 0.06 


0.13 


15 


Wl 


-0.76 ±0.16 


-1.09 ±0.08 


0.33 


31 


W2 


-0.80 ±0.11 


-0.75 ± 0.05 


0.05 


6 


CI 


-0.52 ±0.18 


-0.80 ±0.06 


0.28 


35 



a The reg ions selected for the T-T plots are designated in 
Figure [TBI 

b The errors on qjt are the combination of the 1 — a error 
on the slope of the linear fit and the uncertainty in the flux 
calibration. 



c The median spectral index, ag L , is derived by fitting a 
Gaussi an t o the histogram of the ag values in the region (see 
Figure [TH . The errors are derived from the flux calibration 
uncertainty. 

d Act = — ag^ |. Note that the errors on qjt and ag 

are correlated. 

60% contours have a wider range than on the eastern 
lobe, as at varies from —1.1 to —0.7 with a north-south 
gradient. Southern regions in the tomography maps dis- 
appear into the background for —0.8 < at < —0.7, with 
larger errors on sites close to the edges of the 10% flux 
density contour of the 44 GHz image. Larger spectral 
index variations can be seen in the N-NW sites, with 
—0.9 < at < —0.6. The sites that can be associated 
with spectral index — 1.4 < a t < —1.2 correspond to flux 
density at 44 GHz close to zero and, therefore, are more 
affected by error in the image subtraction. 

The spectral index map, as derived from the direct divi- 
sion of the images at 44 and 18 GHz, is shown in Fig- 
ure QTJ It can be seen that the spectral indices primarily 
vary in the range —0.9 < as < —0.7. In particular, on 
the eastern lobe, within the 90% contours of the peak 
flux density at both 18 and 44 GHz, -0.9 < a s < -0.7, 
while slightly steeper values, —1.0 < as < —0.8, corre- 
spond with the western brightest sites. On both lobes, 
near the northern edge of the 60% contours, the spec- 
tral indices peak at ~ —0.9 on the eastern sites and 
i~ —1.1 on the western ones. Flatter spectral indices 
(—0.6 < as < —0.3) can be located within the central 
and northern regions of the remnant at ~PA 30°. 

It is noted that the regions at the edges of the as map, 
i-e. as < —1.2 (in dark-blue and black) and as > —0.3 
(in light-blue and white) are likely affected by higher 
uncertainty, since these values are associated with lower 
S/N ratios in either the 44 or 18 GHz images. The ef- 
fects of low S/N on the as distribution can be assessed 
from the histogram of the spectral indices over the en- 
tire map. In Figure 1121 the black curve is the Gaussian 
fit over the entire as range, which yields median value 
«s (l = —0.78 and standard deviation a = 0.39. Spectral 
indices as < as —a correspond to regions in the 44 GHz 
image of low flux density (S/N<100), while as > + fl- 
are associated with regions of the 18 GHz image where 
S/N<100. In Figure [T2l the blue curve is a better- fitting 
Gaussian function for as^ — er < as < as^ + o. 



The T-T plots are used to measure the dominant spec- 
tral index within 0'.'5 x 0'.'5 squared regions, as desig- 
nated in Figure Q2] The resultant qtt values are listed 
in Table 21 while in Figure [TJ] the T-T plots for each re- 
gion are shown together with the histograms of the corre- 
sponding as distribution. In seven of the selected areas 
—0.91 < aTT < —0.72, whereas flatter spectral indices 
are derived for the central region CI {oltt ~ —0.52), 
which includes the VLBI position of the SN, and the ad- 
jacent northern region N2 [oltt ~ —0.65). The larger 
differences in the temperature levels within CI and N2 
(see Figures [14"]) . can be seen in the spread of the T-T 
points, which results in the larger error bars attached to 
the linear fit. As shown in in Figurc [T5l the spectral index 
derived for the entire remnant is ctTT m = —0.80 ± 0.05, 
where the flux calibration uncertainty is factored in the 
error (sec §[2]). 

As mentioned in § 15. 1[ among the three methods used 
to assess the spectral index, the T-T plot is the more 
robust, since the spectral index variations resulting from 
the image division and subtraction depend, to different 
extent, on the local intensity of the flux density at the 
two frequencies. However, since the T-T method requires 
boxes larger than the angular resolution of the images, 
the characteristic spectral index of smaller regions, such 
as the features of flatter a t: s at PA ~ 30°, is likely in- 
fluenced by surrounding steeper emission. Otherwise, 
when larger structures are analysed, such as the sites 
of brighter emission on the eastern and western lobes, 
the spectral indices derived from the T-T plots are likely 
more accurate than a local average of at and as- As 
listed in Table 21 over the brightest sites of the remnant 
the discrepancy between axT and regional as^ , which is 
the median value of the as distribution within the T-T 
region, is around 19% over the eastern lobe and peaks to 
~31% in Wl. 

The average spectral index of the western bright regions 
(W1-W2) at ~ —0.78 matches arr m - As discussed ear- 
lier, the steeper as values in Wl are likely affected by the 
higher noise levels in the adjacent northern region. Over 
the eastern lobe, the T-T plots for regions El and E2 
give an average spectral index of —0.86, which is steeper 
than arr m - The steepest value arr ~ —0.91 is associ- 
ated with region El, which covers the brightest sites of 
the remnant. 

From the theory of diffusive shock acceleration (DSA) 
of energetic particles in a uniform magnetic field (see 
iDrurv lfl983l . and references therein), the standard com- 
pression ratio, r, can be derived from the radio syn- 
chrotron spectral index. The axT values given in Tabled 
yield similar compression ratios in regions Wl, W2 and 
E2, specifically with r wi = 2.97±0.43, r W2 = 2.87±0.26 
and TE2 = 2.87 ± 0.14, while a lower compression ratio, 
rgi = 2.65 ± 0.13, results in region El. These results 
suggest a particle spectrum harder over the western re- 
gions than on the brightest eastern sites. For a par- 
ticle density expressed as n(p) oc p~ 7 , where p is the 
scalar momentum, in the scenario of ordinary diffusion, 
n{yv\-wi) °c p~ 2 - 56±0 2e is the spectrum of the western 
side of the SNR, while nisi) cx p- 2 - 82 ±°- 14 j s obtained 
for region El. 
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The lower compression ratio in region El could corre- 
spond with the asymmetric expansion observed in the 
remnant, which appears to be driven by higher east- 
ern v elocities (see § |4j) . As r = pij p\ (| Jones fc Ellison 1 
Il991h , where pi and p\ are the downstream and upstream 
gas densities, respectively, a lower compression ratio over 
the eastern brightest sites could be due to the higher- 
velocity shock drifting the cosmic rays (CR) upstream. 
The CR diffus ion could induce gas heating ups tream 
of the shock (jPtuskin. Zirakas hvili fc Seol 12010(1 and, 
thus, a higher upstream pressure. I f the SNR magneti c 
field is tangled rather than uniform (jDuffy et al. I fl995). 
the softer spectrum could be the result of local sub- 
diffusive particle transport, as the accelerated particles 
get partially trapped in the vicin ity of the shock front by 
struc tures in the magnetic field (jKirk. Duffy fc Gallant I 
1996). In this case, the particle spectrum corresponding 
to region El would be steeper than for diffusive trans- 
port, with rt(ffi) oc p-3-35iEo.i9^ wn ere ry = j DSA (l + 

0.5/r) (IKirk. Duffy fe Gallant lfl996h . 

The softer spectrum in El likely implies that a lower 
fraction of injected particles, i.e. a higher injection ef- 
ficiency, is required for the high emission observed in 
the eastern region. If the stronger eastbound shock 
is coupled with a more efficient particle acceleration, 
the magnetic-field on the eastern lob e mig ht become 
locally amplified (IBell fc Lucek I [20011 lBelTT [20f)l. In 
fact, according to Bell fe Lucek I (|200lh . with efficient 
CR acceleration, (AB/B) 2 ps M, where B and AB 
are the background and fluctuating magnetic fields, re- 
spectively, and A4 is the Mach number of the shock. 
Thus, a higher A4 would likely induce non-linear am- 
plifications of the magnetic field generated in the SNR. 
As the magnetic-field amplification depends on the mag- 
netic field orientation, it can be noted that in asymmet- 
ric bipolar (or bilateral) remnants, such as SNR 1987A, 
regions of higher emission have been linked to cither 
quasi-para llel or quasi-perpendicular magnetic field in- 
clinations (iFulbright fe Reynolds. |[l990UGaensler Ill998t 
lOrlando et al. 112007ft . 

As regards the «rr results in the northern and southern 
parts of the SNR, adjacent small-scale features of both 
steeper and flatter spectra, lead to spectral indices flatter 
than otTT m in the S2 and N2 regions, while the values de- 
rived for regions Nl and SI seem still dominated by the 
steeper indices associated with the bright sites on the 
eastern lobe. We note that, overall, the spectral indices 
over the high-cmissivity regions are steeper than pre- 
dicted by first-order Fermi acceleration at a strong shock. 
This agrees with studies on the spectral distribution of 
young and ra pidly-evolving remnants (e.g. see discussion 
on Cas A in lAnderson et al~~lll99lft . where, since con- 
current high-emissivity and steep spectral indices can- 
not be obtained from first-order acceleration models, 
other mechanisms, such as CR-m ediated shocks and tur- 
bulen t acceleration, are invoked ((Anderson fe Rud nick. 
fl993ft . 

In both at and as distributions, gradual transitions 
from steeper to flatter spectral indices (i. e. —0.8 < 
o^t.s ^ —0.5) can be seen on the edges of the rem- 
nant and around, or in proximity, of the brightest 
sites. Spectral index gradients in older and larger 



SNRs have been explained as d ue to multiple overlap- 
ping spectral index structure s (|DeLanev et al. I 120021 : 
iTam. Roberts fc Kaspi. I I2002T) . Accordin g to hydrody- 
namic simulations coupled to DSA by lEllison et al. I 
(|2007ft of a spherical SNR, a spectral gradient might 
reflect a spectral index structure over multiple con- 
centric shells of shocked material within the forward 
and the reverse shocks, with a CR population vary- 
ing between the concentric shells, such that the spec- 
tral index results flatter on the outer shell than on 
the inner ones. The scenario of sub-diffusive trans- 
port of particles at the shock front, also leads to 
small-scale spectral gradients in the downstream plasma, 
which would correspond to the transition between sub- 
diffusi ve and diffusive behaviour of the particle propaga- 
tion (jKirk. Duffy fc Gallant lll996l ). 

All methods used for spectral analysis show two inner 
regions of flatter spectral indices, where — 0.6 < axr ^5 
—0.5 and —0.5 < Ot,s —0.3. The more central fea- 
ture, which overlaps with the SN site (see Figure ITTj) . has 

a t « -0.5, -0.6 < a s < -0.5 and a TT 0.52 ± 0.18. 

The second feature, located further north, can be asso- 
ciated with at ~ —0.4 and —0.4 < as < —0.3, while 
the corresponding aTT, between regions Nl and N2, is 
affected by the nearby steeper spectra. 

We note that a central feature of flatter spectrum has 
been seen in other two-frequency spectral index image s 
of the SNR (jPotter et al. \mM . iLakicevic et al. \\20m . 
In Figure[TT]the spectral index image is also overlaid with 
the contours of the HST observations in 20110 (in black), 
and with the 10% level of the peak flux density of the 
high-resolution image at 44 GHz (in yellow). It can be 
seen that the central feature of flatter electron spectrum, 
not only coincides with the inner region of fainter emis- 
sion that is visible in the radio at higher resolution, but 
also partially overlaps with the western side of the optical 
ejecta. On the contrary, the northern feature, which ap- 
pears to connect with the central feature along PA ~30° , 
approximately in alignment with the main direction of 
the optical ejecta, does not seem to have been noted be- 
fore. This feature exhibits at,s values flatter than the 
central region, and, since it appears to be smaller than 
the O^ x 0'.'5 boxes used for the T-T plots, cannot be 
associated with a local «tt value. Nevertheless, since it 
corresponds to a site characterised by high flux density 
both at 18 and 44 GHz, the related at and as values are 
unlikely affected by significant error. 

While pulsar wind nebulae (P WNe) in the radio ban d 
are identified by -0.3 < a < (jGaensler fc Slane II2006T ). 
the acceleration origin in the above discussed features of 
—0.5 < a < —0.3, could be explained in terms of injec- 
tion by a pulsar. Further observations at high frequencies 
and monitoring of any morphology changes is required to 
clarify the nature of these features. 

6. EMISSION AT OTHER WAVELENGTHS 

As sh own in early models of Type II SNRs ((Chevalier I 
11982ft . the SN-CSM interaction in SNR 1987A has gen- 

2 STScI-2011-21, NASA, ESA, &; Challis P. 

(Harvard-Smithsonian Center for Astrophysi cs) 

http : //hubblesite . org/newscenter/archive/releases/2011/21/i 
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erated a main double-shock structure, which includes, 
from the outside inwards, the forward shock and the re- 
verse shock. Between the forward and reverse shocks, 
reflected shocks should be present due to the forward 
blast wave encountering the high-density CSM in the 
ER (jBorkowski. Blondin fc McCrav I [l997h . While the 
structure of the shock can be probed by comparing the 
remnant morphology in the radio and X-ray wavelengths, 
the distribution of the progenitor dense material can be 
investigated via Ha images. 

The radio emission mainly originates at the forward 
shock as this collides with the dense CSM associated with 
the ER, thus creating a discontinuity i n the magnetic 
field where particles are accelerated (IZanardo et al. I 
2010). The forward shock can be identified with the 
sharp outer edge of the SNR shell, which propagates 
into the CSM of spatially varying density. A similar sce- 
nario happens at the reverse shock, or inner edge of the 
SNR envelope, which is normally dominated by the X- 
ray emission. Between the reverse shock and the ER, 
the reflected shocks lead to the higher temperatures of 
the shocked gas and a flatter density profile. Up to day 
^4000, the X-ray emission has likely been generated by 
the interactions with th e low-density Hli region locate d 
on the inside of the ER ( [Chevalier fe Dwarkadas 1 1995). 
whereas, in the last ^5000 days, it has primarily orig- 
inated from the i nteraction with the dense inner ring 
(jPark et al. II2011D . 

In Figure [16j the 4 4 GHz image is com pared to recent 
X-ray obs e rvatio ns (|Helder et al. I [2012). As noted by 
INg et al. I (HSU, the X -ray emission exhibits an east- 
west asymmetry less marked than what seen in the ra- 
dio, while the north-south asymmetry is significant. In 
particular, the brightest X-ray site, which appears in the 
NE quadrant between PA 30° and 60°, overlaps with the 
brighter regions on the eastern lobe of the 44 GHz im- 
age, but peaks slightly northwards of the eastern radio 
peak. In the western lobe, the X-ray emission peaks at 
PA -270° and between PA 300° and 330°. While the 
radio image is also reaching the western peak of bright- 
ness at PA -270°, the region between PA 300° and 330° 
corresponds to one of the fainter sites (see § 13.11) . The 
offset between the X-ray and radio peaks is demonstrated 
in the polar projection of the two images shown in Fig- 
ure nu 

The 44 GHz map and the HST image from recent obser- 
vations are superimposed in Figure 1171 It can be seen 
that the morphology of the radio emission matches the 
ring optical features. As shown in Figure 111! on both 
the east and west sides of the optical ring, the eastern- 
most and westernmost hot spots appear to overlap with 
the bright regions on the eastern and western lobes of the 
radio images. The fainter regions of the ring visible at 44 
GHz, i. e. located in the NW and SW quadrants, seem to 
coincide with sites where clusters of Ha hot spots appear 
somewhat disjointed. This is especially pronounced on 
the NW quadrant, between PA 300° and 330°, where the 
forward shock might have overtaken the ER and reached 
less dense CSM regions, while the inner surface of the 
ring is still hot enough to emit X-rays. 

As regards the remnant asymmetry, contrary to the ra- 
dio images, the optical emission now appears markedly 



brighter on the western lobe. The difference in the asym- 
metry direction between the nonthermal radiation and 
the Ha emission is more clear in the RGB overlay shown 
in Figure [THJ It can be noted that the eastern lobe of 
the ring is dominated by the radio (red) and X-ray (blue) 
emissions, which turn into shades of violet in the overlay, 
especially in the NE quadrant where the radio and X-ray 
brighter sites overlap. On the western lobe, the fainter 
radio emission is now overshadowed by the Ha emission 
(green). This appears to validate the hypothesis that 
the remnant asymmetry is likely due to an asymmetric 
explosion of the progenitor rather than to an asymmet- 
ric distribution of the CSM. We may therefore expect 
the radio and X-ray emission to follow suit, with the 
east-west asymmetry gradually reversing. The possibil- 
ity of an asymmetr i c init ial explosion was suggested by 
iChevalier fc Sokerl (|1989f) to explain the asymmetric ex- 
pansion of the remnant, whereas hydrodynamic calcu- 
lations of the remnant evolution from the asymmetric 
explosion of a ~20 M© merger of binary systems, have 
been proposed to explain the BSG progenitor and to link 
the BSG-RSG wind interaction to the triple- ring neb- 
ula structure ( Podsiadlowski. Morris fc Ivanova. 1 120071 : 
iMorris fc Podsiadlowski 1120071 12009| )~ 

In Figure 1171 the radio emission from the inner re- 
gion of the remnant as seen at 44 GHz can be com- 
pared to the optical emission from the ejecta. It is noted 
that the ejecta is blue an d red-shifted along the line- 
of-sight (jKjaer et al. I [2010) and, like the emission from 
the ER, the optical ejecta is characterised by an east- 
west asymmetric mor phology, whose origin is yet unclear 
(|Larsson et al. Il2011f) . The western side of the ejecta ap- 
pears to overlap with the eastern part of the central fea- 
ture visible at 44 GHz and, as mentioned in § 15. 2\ extends 
over a region of the SNR that can be associated with flat- 
ter spectral indices. Moreover, the 'hole' in the ejecta, 
which is close to the SN location (see Fig ure ITTj) and has 
becom e more pronounced in recent years (|Larsson et al. I 
1201 If ) , also overlaps with the central emission detected at 
44 GHz. 

7. FREE-FREE ABSORBED COMPACT SOURCE 

iMatsuura et al. I (|2011| ) have estimated that the largest 
dust component, with mass 0.4 — 0.7tfo, sits in t h e cen- 
tral region of the SNR. However Lakiccv ic et al. I (|2012f ) 
have concluded that, for observations up to 94 GHz, the 
dust does not contribute to the radio emission. The cen- 
tral feature detected at 44 GHz might therefore originate 
from a compact source or PWN located in the vicinity 
of the SN ejecta. 

To determine whether non-thermal radio emission from 
a pulsar or a PWN located in the inner regions of the 
SNR would be detectable, we estimate the fraction of 
this emission that would be absorbed by the ionised gas 
likely located within the ER. Assuming the ionised mass 
in the ejecta, M e j, to be distributed within a sphere 
of radius R, the free- free optical depth, r//, is pro- 
portional to the path length through the free-free ab- 
sorbing spherical region along the line-of-sight (los). In 
convenient astronomical un its, tjj can be estimated as 
(jRvbicki fc Lightman Ill979l) 

r ff ^T T- 3/2 ^- 2 gJjEM 
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where T is the temperature of the plasma in K, v is the 
frequency of the flux density in GHz, gJJ = gjj(v^ T) 
is the velocity-averaged Gaunt factor appropriate to the 
observing frequency, and the emission measure is defined 
as 

EM = I n e rii dl 

J\os 

where n e and rii are the electron and ion densities, re- 
spectively, expressed in cm -3 , and I is the path length, 
expressed in pc, along the line-of-sight. As Hn regions 
typically have electron temperatures of order 10 4 K, in- 
dicated by T4, and, in the radio regime, gJj(h , ,T) w 
5.99T 4 015 i/- 01 (jBrown 1119871 ) . the free-free opacity be- 
comes 

TffK 3.28 x 1Q- 7 T^ 135 iy- 21 nlR 

assuming n e — m, which is approximated by 
rie(M e i R ) sa 3 M ej (m p 4ir R 3 )' 1 , and Jj o dl « R. 
Figure [19] shows the variation of 77 / as a°!unction of the 
ionised mass for M e j < 10 Mq and volume size R < 0.1 
pc. The contour levels are plotted for < 77/ < 2. 
According to this simplified model, an ionised mass in 
excess of 5 Mq within a radius of 0'.'2 (approximately cor- 
responding to the size of the densest part of the ejecta 
visible with HST), is required to produce a free- free op- 
tical depth of unity at 44 GHz. We note that the flux 
density from the emitting plasma, within a ~ 0'.'2 solid 
angle, given by SM ej — (2fcsT/A 2 )ri, is negligible. In 
this scenario, the radiation observed through the region 
of ionised plasma would be S b s ~ 0.63 S , 44 C , where S44 C 
is the actual flux density at 44 GHz associated with the 
compact source in the background. Therefore, since the 
flux density measured in the inner region of the 44 GHz 
image is S b s ~ 1.4 mJy (see § 13.11) . the actual radio 
emission originating from a compact source inside the 
ER could be as high as ~ 2.2 mJy. Using the spectral 
index derived in Figure [6j this emission would lead to 
Sq4 C ~ 1.3 mJy at 94 GHz, which is consisten t with the 
upper limit of 1 mJy set by lLakicevic et al. I (|2012j ) at 
that frequency, for any discrete radio source in the cen- 
tral region of the remnant. 

In the case of a more r ealistic, steeper density p rofile, 
such as that proposed bv lChevalier fc Soker I ()1989l ) with 
n e oc R~ 9 6 , the emission from a PWN would be more 
heavily absorbed and could be detected only if there are 
holes in the ejecta. As it has been noted from HST 
images in the R-band and B-band , the ejecta is likely 
to have holes (jLarsson et al. l[20lTh ( see Figures [T7] and 
|T5J) , thus the radiation emitted by a PWN or compact 
source might escape the dense ionised material and be- 
come visible at lower frequencies. On the other hand, 
with a more compact density distribution, part of M e j 
may have fallen back onto the neutron star that formed 
in the explosion, and, subsequently, the neutron star and 
any accreted ejecta might have already collapsed into a 
black hole. 

8. CONCLUSIONS 

We have presented the first image of SNR 1987A at 
44 GHz, as derived from ATCA observations performed 
in 2011. To investigate the spectral index distribution 
across the remnant, this new image has been analysed 



in conjunction with that derived from contemporaneous 
observations at 18 GHz. The emission morphology has 
been also compared to contemporaneous optical and X- 
ray observations. A summary of our findings is as fol- 
lows: 

• Consistent with previous radio observations, the 
44 GHz image shows a marked asymmetry in the 
emission distribution. More specifically, the east- 
west asymmetry ratio is ~ 1.5 from the ratio of 
the brightness peaks in the radial profiles at PA 
^90°, is ~ 1.6 from the integrated flux densities 
over the eastern and western halves of the image. 
These values are higher than the ~1.4 ratio derived 
for the new 18 GHz image and the ratio previously 
measured with images at lower frequencies. 

• The comparison between the new images at both 
18 and 44 GHz with corresponding observations 
performed in earlier epochs, specifically the 2003 
observations at 17 and 19 GHz and the 2008 obser- 
vations at 36 GHz, highlights an asymmetric ex- 
pansion of the remnant, with expansion velocities 
on the eastern lobe significantly higher than what 
measured on the western lobe. 

• The 18-44 GHz spectral index distribution is mea- 
sured at an angular resolution of 0'/4. The spec- 
tral indices in SNR 1987A primarily range between 
—1.1 and —0.3, with a mean of —0.8. Spectral in- 
dices associated with the brightest sites over the 
eastern lobe are steeper than the mean value. The 
steeper spectrum on the eastern lobe implies com- 
pression ratios slightly lower than on the western 
bright sites, and could be correlated with the higher 
expansion rate measured on the eastern side of the 
remnant. Two regions of flatter spectral indices are 
identified, one approximately located in the centre 
of the SNR and the other located further north. 
These two features lie at PA ~ 30°. 

• There is a strong correspondence between major 
features of the emission at 44 GHz, and the ar- 
rangement of the hot spots shown in the Ha emis- 
sion. The direction of the east-west asymmetry of 
the X-ray and radio emission, is opposite to that of 
the Ha emission. This fact supports the hypothesis 
that the remnant asymmetric morphology might be 
due to an asymmetric explosion, rather than to an 
asymmetric distribution of the CSM. 

• At 44 GHz, a central feature of fainter emission 
appears to extend over the SN site, and to overlap 
with the western side of the ejecta as seen by HST. 
This feature corresponds to a region of flatter spec- 
tral indices in the 18-44 GHz spectral map, which 
could indicate the presence of a compact source 
or a PWN. The origin of this emission is unclear. 
However, simple free-free absorption models sug- 
gest that the radiation emitted by a compact source 
inside the equatorial ring may now be detectable at 
high frequencies, or at lower frequencies if there are 
holes in the ionised component of the ejecta. Fu- 
ture high-resolution observations, both at lower fre- 
quencies with VLBI and at higher frequencies with 
ATCA and the Atacama Large sub-Millimeter Ar- 
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ray, will be crucial to further investigate the nature 
of this emission. 
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5 h 35 ra 28.2 s 38.1 s 28.0 s 27.9 s 27.8 s 

Right Ascension (J2000) 



Fig. 1. — Top: Diffraction-limited Stokes-I continuum image of SNR 1987A at 44 GHz made by combining observations performed with 
the ATCA on 2011 January 24 and 2011 November 16. The beam size is 0'.'35 X 0'.'23 as plotted in the lower right corner. The off-source 
rms. noise is ~0.02 mjy beam -1 . Bottom: Slightly super-resolved 44 GHz image obtained using a 0'.'25 restoring circular beam (lower 
right corner). 
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5 h 35 m 28.3 s 28.2 s 28.1 s 28.0 s 27.9 s 27.8 s 27.7 s 

Right Ascension (J2000) 




Right Ascension (J2000) 

Fig. 2. — Top: Diffraction-limited Stokes-1 continuum image of SNR 1987A at 18 GHz from observations performed on 2011 January 26 
using the ATCA. The beam size is 0'.'63 X 0'/47 as plotted in the lower left corner. The off-source rms. noise is ~90 /^Jy beam -1 . Bottom: 
Super-resolved 18 GHz image obtained using a 0'.'25 restoring circular beam (lower right corner). 
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Fig. 3. — Radial slices through the diffraction-limited 44 GHz image at 6 position angles. Black dots on vertical axes indicat e the 
position of the zero for each slice. The offset is radial from the VLBI position of SN 1987A as determined by Reynolds ct al. ( 1993) [RA 
05 h 35 m 27?968, Dec -69° 16' 11'.' 09 (J2000)], and is positive toward east and/or north. 
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5 h 35 m 28.2 s 28.1 s 28.0 ! 27.9 s 27.8 s 

Right Ascension (J2000) 

Fig . 4. — Comparison between the 44 GHz image and that derived from ATCA observations at 36 GHz performed in 2008 ( Potte r~et al. I 
12009 ). The 4 4 GHz image is rest ored with the 0'.'33 X 0'.'24 beam that characterises the diffraction-limited image at 36 GHz, which is set 
at PA -1.3° ijPotter et al. Il2059l . TableEJ. The contours at levels 25% - 85% of the peak flux density, in step of 10%, are in blue for the 
36 GHz image and in orange for the 44 GHz image. To highlight the morphological differences of the brightest sites in the two images, the 
44 GHz contours at 45%, 65 % and 85% of the peak flux density, are filled in yellow, orange and brown. The cross-hair symbol marks the 
VLBI position of SN 1987A HRevnoids et al. IITS95T) . 
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Right Ascension (JKOOO) 

Fig. 5. — Comparison of the diffraction-limited 44 GHz image (top), as derived from ATC A observations pe rformed in 2011 January and 
November, with the image generated with the Fourier-domain model (centre) developed by |Ng et al. I (2008) to fit the SNR observations 
at 9 GHz. The model used in the comparison fits the 9 GHz observations performed in 2011 April 22 (Ng et al., in preparation), and has 
a radius of 0'.'917. To match the total flux of the observations at 44 GHz, the model flux density has been scaled using a spectral index 
a = —0.74 (see Figure[BJ. Bottom: Dirty map of the residual visibilities obtained by subtracting the model from the 44 GHz observations. 
No deconvolution has been applied to the residual. 
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Fig. 6. — Spectral index of SNR 1987A as determined from ATCA observations performed in 2011. In detail, the spectral index is 
determined from the power-law fit to the flux densities measured at five frequencies: 1.4 GHz (Zanardo et al., in preparation) (solid circle), 
9 GHz (Ng et al., in preparation) (hollow triangle), 18 and 44 GHz (this paper) (hollow circle), and 94 GHz (Lakicevic et al. 2012) (solid 
diamond). The fit yields a = -0.74 ± 0.01. 
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Fig. 7. — Radial slices through the diffraction-limited 18 GHz image at 6 position angles. Black dot s on vertical axes indi cate the position 
of the zero for each slice. The offset is radial from the VLBI position of SN 1987A as determined by Reynolds ct al. (1995), and is positive 
toward east and/or north. 
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5 n 35 m 28.2 



2B.1 5 28.0 s 27.9 s 

Right Ascension (J2000) 



Fig. 8. — Comparison between the new 18 GHz image and that derived from ATCA observations at 18 GHz performed in 2003 
(Manchest er et al. II2005T ). Both images are super-resolved with a 0'.'25 circular beam. The contours at levels 25% — 85% of the peak 
flux density, in step of 10%, are in blue for the 2003 image and in orange for the new image. To highlight the morphological differences of 
the brightest sites in the two images and the significant expansion of the remnant since 2003, the contours of the new image at 45%, 65% 
and 85% of the peak fl ux density, are filled in yellow, light and dark orange. The cross-hair symbol marks the VLBI position of SN 1987A 
(Reynolds et al. 1995). 
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Fig. 9. — Remnan t expansion from images at 18 GHz on day 6003 (Manchester ct al. 2005) and day 8738, at 36 GHz on day 7815 
UPotter ct al. 2009) and 44 G Hz on day 8886 (this paper). The expansion is estimated from the distance between the VLBI position of 
the SN ( Reynolds ct al. 1995) and the position of the peak of brightness for the eastern and western lobes, as derived from the PA 90° 
emission profiles centred on the SN coordinates. The blue symbols are associated with the brightness peaks on the eastern lobes, while the 
red ones represent the distance between the SN position and the brightness maximum in the western lobe. Mean values derived for each 
pair of measurements are represented by the black symbols. The blue line is the linear fit of the measurements at 18 GHz from 2003 to 
2011, while the purple line only fits the 36 and 44 GHz within 2008 and 2011. For data between 20 08 and 2011, the derived expansion 
velocity is 3900 ± 300 km s _1 . As a comparison with the expansion rate derived by Ng ct al. (2008), the dashed black line is the slope 
of the linear fit of radius measurements extracted from the 1992-2008 super-resolved images at 9 GHz fitted with a torus model, the error 
associated with the fit is represented by the grey-filled area. 
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Fig. 10. — Tomographic spectral images of the remnant of SN 1987A, where the spectral index is defined as S oc u a . The spectral images 
are derived from the 18 GHz image from observations performed on 2011 January 26, and the 44 GHz image from observations between 
in 2011 January and November. Both images are centred on the VLBI position determined by Reynolds et al. (1995), and restored with a 
circular 0'.'4 beam (top panel, lower left corner). The grey scale shows difference images for —1.4 < at < —0.3, in steps of 0.1. Blue/white 
contours correspond to the 44 GHz data convolved with a C^'4 circular beam at levels 60%— 90% of the peak flux density, while the related 
contours corresponding to the 18 GHz image are in red. In the top left panel, axes at six position angles, centred on the SN position, are 
also indicated. All images are set on the same intensity range. 
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Fig. 11. — Top: Map of the 18 — 44 GHz spectral index distribution. Both the 18 and 44 GHz Stokcs-I images have been restored with a 
0'.'4 circular beam, and regriddcd at a pixel scale of 3 mas. Image regions below the 44 GHz rms noise level were masked. Contours of the 
60% and 90% flux density levels for the 18 GHz {yellow) and 44 GHz (blue) images, identify the emission peaks on each lobe. Bottom: To 
locate the inner regions of the spectral index map with respect to the emission measured at 44 GHz, the outline of the diffraction-limited 
image is shown via its contour at 10% of the peak flux density (yellow). To locate the density distribution of the CSM, the spectral index 
map is also overlaid with the contours of the H ST image derived from observations performed in 20110, specifically at the 2%, 7%, 12% 
and 20% emission levels (black). As indicated by Larsson ct al. (2011), the innermost ffi>T contour of the ejecta (gray fill), east of the SN 
position (black star), corresponds to the location of an emission drop, or 'hole', in the optical image. 
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-1.5 -1.0 -0.5 0.0 0.5 

Spectral Index 



Fig. 12. — Histogram of the spectral index values, as, derived from direct division of the images at 44 and 18 GHz (see Figures [TT} . 
The vertical axis is the normalised number of occurrences. The black line is the Gaussian fit of the whole histogram, with median value 
ctg^ = —0.78 and standard deviation a = 0.39. The blue line is the Gaussian fit for ag^ — a < as < ag^ +a. Spectral indices ag < ag^ —a 
correspond to regions in the 44 GHz image of low flux density (S/N<100), while spectral indices a$ > ag^ + a are associated with regions 
of the 18 GHz image where S/N<100. 



High-resolution radio observations of SNR 1987A at high frequencies 



23 




28.10s 28.00s 27.90s 

RA (J2000) 



27.80s 5h35m27.70s 



Fig. 13. — Two frequency (18-44 GHz) spectral index inten sity map with indication of nine 0'.'5 X 0''5 regions (black squares) used for 
the temperature-temperature (T-T) plots presented in Figure [14] For each region, the spectral index values, «tt, deri ved from the T-T 
plots, are noted in yellow, while the median values derived from the spectral map, aj , are in white font (see Figure ITU . The contours of 
the 60% and 90% levels of the peak flux density of the 44 and 18 GHz images, resolved with a 0''4 circular beam, are shown in blue and 
red, respectively. 
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Fig. 14. — Two frequency (18-44 GHz) temperature-temperature (T-T) plots for nine remnant regions, designated as in Figure [l3l For 
each region, the histogram of the distribution of the spectral index values, ag, is also shown, as obtained from the spectral index intensity 
map. In each histogram figure, the solid black circle indicates the related median spectral index and the dashed black line corresponds to 
the standard deviation, as derived from the Gaussian fit of the histogram. The solid blue circle indicates the spectral index derived from 
the linear-regression fit of each T-T plot, aTTi ancl the blue bar indicates the 1 — tr error. The values of chj-t an d a S li associated with 
each region are listed in Table [4] 
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Fig. 15. — Temperature-temperature (T-T) plot of the whole SNR at 44 GHz and 18 GHz (black squares). The best fitting straight 
line corresponds to a = —0.80 ± 0.01(±0.05) (S = v a ), where the two errors represent random and absolute flux density scale errors, 
respectively. 
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Fig . 16. — Left: Overlay of the Chandra 0.5-8 keV image (contours) derived from observations performed in September 2011 (Helde Fet al. I 
120121) . on the 44 GHz diffraction-limited image from 2011 ATCA observations. The X-ray contours are at levels 20%— 80% of the emission, 
in steps of 20%. Right: We show the 44 GHz diffraction-limited image converted to polar coordinates to visualise asymmetries in the 
2-dimensional radial distribution. In this image, the position angle and the projected radial distance from the geometrical centre of the 
remnant, are the new coordinates. The emission intensity is conserved in the conversion. 
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Fig. 18. — RGB overlay of the optical image from 2011 HST observation^] (green), the X-ray image from Chandra observations performed 
in 2011 (Hcldcr et al. 20l2) (blue) and the 44 GHz image from ATCA observations in 2011 January and November (red). 
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Fig. 19. — Variation of the optical depth, tjj, of an idealised Hn region in SNR 1987A, approximately located at the ejecta site. The 
ionised mass in the region is assumed to vary in the range < M e j < IOMq , and to be uniformly distributed within a spherical volume of 
size < R < 0'.'4, which corresponds to a path length up to 0.1 pc along the line-of-sight. The yellow colour indicates low free-free opacity; 
the brown colour indicates high opacity. 



